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Abstract We make some remarks on reconnection in plasmas and want 
to present some calculations related to the problem of finding velocity fields 
which conserve magnetic flux or at least magnetic field lines. Hereby we 
start from views and definitions of ideal and non-ideal flows on one hand, 
and of reconnective and non-reconnective plasma dynamics on the other 
hand. Our considerations give additional insights into the discussion on 
violations of the frozen-in field concept which started recently with the 
papers by Baranov & Fahr (2003a;2003b). We find a correlation between 
the nonidealness which is given by a generalized form of the Ohm's law and 
a general transporting velocity, which is field line conserving. 



1 Introduction 

For many applications in astrophysics it is interesting to ask by which veloc- 
ity the magnetic flux is transported. This problem of MHD subject was al- 
ready analysed in several articles, e.g. Newcomb (1958), Vasyliunas (1972), 
Schindler & Hesse (1988), Hesse & Schindler (1988), Hornig & Schindler 
(1996). Here we shall make use of the basic ideas of these articles and apply 
these to the corresponding more specific problem in hcliospheric interface 
physics. If no appropriate velocity field for the magnetic flux transport can 
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be found in a plasma region then one can identify the occurence of magnetic 
reconnection. The existence of a flux-transporting velocity field is connected 
with the 'ideal' constraint, that V x R = 0, where R = E + v x B de- 
notes the nonidealness, E, v and B being the electric field, the plasma 
bulk velocity, and the magnetic field, respectively. The above requirement 
should be fulfilled everywhere in the regarded domain, with the possible 
exception of some localized regions of non-vanishing resistivity. Inside the 
resistive domain the plasma velocity needs not to be flux-conserving, but, 
under certain circumstances, a velocity field different from v could perhaps 
be found which does the freezing-in of B. Such a new velocity field therefore, 
should be continuous over the whole domain and converge to the normal 
flux conserving plasma velocity outside the resistive domain. 

If it is not possible to find such a new velocity field, w (where w = flux 
velocity), then magnetic reconnection is taking place since the flux trans- 
porting velocity field is not continuous at the border between the ideal and 
the nonideal region. This is shown schematically in Fig.^ where the flux ve- 
locity, w, within the nonideal region is different from the plasma velocity, v. 
Outside the nonideal region, the flux velocity, w, equals the plasma velocity, 
v, as shown in Fig.[3 Therefore, the nonideal region seems to have discon- 
nected flux tubes. This means that every magnetic reconnection process is 
a nonideal process, but not vice versa. 

Here we want to answer the question, how the MHD plasma velocity field 
can be redefined, in such a way, that it is again the virtue of a flux- or at least 
a line-conserving flow. One actual reason for that discussion is the recent 
dispute between Baranov & Fahr (2003a; 2003b) and Florinski & Zank 
(2003). These authors discussed the meaning of a generalized Ohm's law, 
where the nonideal part describes the interaction of different particle species 
in a partially ionized plasma. Baranov and Fahr claimed, that the magnetic 
field, under specific conditions valid in the heliospheric interface region, 
may neither be frozen into the mass-weighted plasma bulk flow, nor in the 
ion velocity. We can show very easily, that with the assumptions and the 
form of Ohm's law the authors have found (Baranov & Fahr, 2003a), a flux 
conserving form of Ohm's law could in principle be constructed. However, 
the needed flow velocity is not necessarily identical with one of the species 
bulk velocities, as proposed by Florinski & Zank (2003), but critisized by 
Baranov & Fahr (2003b). Florinski & Zank claim, that the magnetic flux 
is frozen into the ion velocity within and in the vicinity of the heliosphcre. 
Baranov and Fahr have shown, that in fact this depends very sensitively on 
the solution of the whole set of the multifluid equations, and especially on 
the magnetic field structure. 

We want to show how flux or line conserving velocities could be calcu- 
lated for a general form of an Ohm's law. In addition we discuss how the 
technical procedure is practised and how the included considerations can 
principally be used to solve this problem for a partially ionized plasma like 
that of the heliospheric interface. 
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Fig. 1 The part of a field line (magnetic flux tube) that is penetrating a nonideal 
region (shown as the grey shaded sphere) experiences a different velocity than the 
rest of the field line that is still outside. The figure is taken from Priest et al. 
(2003). 




Fig. 2 Field lines passing through a nonideal region can disconnect. Outside 
the nonideal region they move with the same smoth velocity field, w is the flux 
transporting velocity, v is a species bulk flow velocity or plasma velocity and tilde 
means vector. The figure is taken from Priest et al. (2003). 

2 Derivation of flux conserving velocity fields 

There arc two basic concepts of ideal, respectively non ideal processes in 
MHD: that of line conservation (violation), which marks a breakdown of 
magnetic line connectivity, and that of magnetic flux conservation. If a 
highly localized breakdown or violation process takes place in an ideal 
plasma environment, this is called magnetic reconnection. However, even 
in a nonideal plasma it is possible to get solutions of the nonideal MHD 
equations, conserving magnetic flux. For this to happen it is only necessary 
that 

VxR = (1) 

for any arbitrary closed fluid line (see Priest & Forbes, 2000). The nonide- 
alncss R is given by 

R=E+vxB (2) 

An explicit discussion of nonideal terms can be found, e.g. in Schluter (1958). 
The second concept of a frozen-in magnetic field is not as strict as the upper 
one and only requires 

B x (V x R) = (3) 



which is only fulfilled, if 

V x R = AB 



(4) 
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where A is a function of location and time in general and B • VA = 0. If 
therefore 

B x V x (E + v x B) / (5) 

is valid, magnetic field line connection is breaking down and magnetic re- 
connection is taking place. 

In the stationary case (V x E = 0) the criterion can be written as 

B x V x (v x B) ^ (6) 

For this reason we want to find a flux conserving velocity field w and a 
function I with 



E + vxB = R <a> E|wxB = VI (7) 
Eq. J2J can be written as 

(w - v) x B = VI - R (8) 
Therefore we can formulate 

Theorem 1 : If Ohm 's law is given by E + v x B = R, where v is the 
plasma velocity, then the magnetic flux is frozen-in with respect to velocity 
fields 



B x (VI - R) 
IBP 



M B (9) 



or, equivalently 



(E - VI) x B 
w=^ — J- + /2B (10) 



where fi, i.e. fi is any arbitrary function in space and X is the solution of 
the partial differential equation 

B • VI = R • B (11) 

with appropriate boundary conditions. This procedure is reasonable, although 
it can be seen in the above derivation that even in the ideal case there exists 
an infinite number of alternative velocity fields in all of which the magnetic 
field could be frozen in. This additional velocity component is directed in 
the 'binormaP direction of the magnetic field, writing 8v | B x VI. This 
follows due to B • VI = 0, so I is an integral of B (in the ideal case). That 
field line motion is not unique, not even in the case of an ideal plasma, was 
also discussed, e.g. in Vasyliunas (1972). The bulk velocity field or plasma 
velocity can therefore be regarded as a minimal flux preserving velocity (see 
Vasyliunas, 1972). Therefore the function I in equation should be cho- 
sen carefully (also the term parallel to the magnetic field fiB). With respect 
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to physical interpretation w should be a velocity field which has a reason- 
able physical meaning. Under certain circumstances, the magnetic flux in 
a partially ionized plasma (see the discusssion in Baranov & Fahr (2003a; 
2003b)) nearly is frozen into the ion velocity. On the other hand, in a nearly 
completely ionized plasma electrostatic turbulences could lead to a strong 
localization of the nonidcalncss (or resistivity in certain cases), so that the 
flux transporting velocity should at least smoothly converge into the normal 
known plasma bulk velocity outside or to say far away from the localized 
nonideal region, which is typical for astrophysical plasmas. This can be 
guaranteed by the boundary condition 

X = constant Vx^D\Dr and especially on DDr (12) 

where Dr is the nonideal domain. If, however, the magnetic field is nowhere 
frozen in the bulk fluid motion, then one has to drop condition (|12H , and the 
plasma is everywhere frozen in the velocity field @ • Then the nonidealness 
is not localized, and one can hardly speak of a localized nonideal instability 
or a localized reconnection process. 

Let us now consider a situation in which R vanishes outside a certain 
domain Dr d D (or goes faster to zero or to a certain limit). Then we can 
take equation Qlljl- and use the identities 

E B = R B (13) 
V • (IB) = B VI (14) 

to write down 



D 



B-VXdV = J R BdV (15) 
b 

^ J V-{XB)dV = J E-BdV (16) 

D D 

J XBdS = J E-BdV (17) 

dD D R 

A similar discussion was done by Priest et al. (2003). These authors empha- 
sized the importance of the component of the electric field Em being aligned 
with the magnetic field. Here we take a look to this problem from a different 
point of view, which is also connected with E|| , but emphasizing the role of 
the Dirichlet boundary condition for X. From this we can see that taking 
the boundary condition (fT2"|l , there is only a possibility to solve Eq. (fTTJ) , if 
and only if the right hand side integral over the domain Dr vanishes identi- 
cally. This would imply, that the different field aligned parts of the electric 
fields have to cancel out 'statistically' , e.g. as a result of a strongly spatially 
fluctuating electric field component due to turbulence with no preference in 
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direction with respect to the magnetic field. If there is no preferred direc- 
tion inside the resistive domain there is at least one surface with vanishing 
E ■ B = 0, where the Lorentz invariant changes its sign. If, in contrast, this 
integral on the right hand side of Eq. (|17fl does not vanish and has a non 
negligible value, there is no possibility to find a velocity field which exists 
everywhere in the whole domain and which is continuous across the border 
of the resistive region Dr. Field lines crossing both, the resistive and the 
ideal region, are convected with a velocity field, which is not identical with 
the plasma velocity, even outside Dr, and this implies serious breakdown of 
magnetic flux conservation. The reason for this is, that field lines being out- 
side and crossing Dr have different flux velocities compared to field lines, 
which do not cross Dr and being completely outside Dr. That is, what can 
be called reconnection. 



3 Derivation of line conserving velocity fields 

Theorem 2 : If E + v x B = R is generalized Ohm's law, then for 
ordered, non-ergodic fields, it is always possible for each solution of the 
above equation to find a velocity field w and a function A so that 

BxVx(E + wxB) = BxVxY = (18) 
V x Y = AB (19) 



and 

B x (Y - R) 



where Y has to fullfill the equations 



AiB (20) 



B (R Y) = (21) 
VxY = AB => B-VA = 0. (22) 

The existence of fields Y, fulfilling Eq. \2ty. together with the existence 
of a corresponding field line constant A imply magnetic field line conserva- 
tion (see e.g. Vasyliunas, 1972). Line conservation means, that two fluid 
elements are always connected by one field line during the convective evolu- 
tion of the electromagnetic field and the velocity field. 

Proof: 

We are searching for 

E + vxB = R <=>■ E + wxB = Y (23) 
=> (w-v)xB = Y- R (24) 

If one inverts the left hand side of Eq. I|24|l we get Eq. i|2U|) and as neces- 
sary additional conditions for Y Eq. I|21l) with line conservation condition 
Eq. 123 . 
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Do these fields A, Y exist? If we use Euler-potentials 
Y = Y a Va + Y^VP + Y s Vs 
and take the curl of Eq. (|25|> , we get 



(25) 



V x Y = 



dY fi q Y <* 

w 



da 



ds 



(Vq x V/3) + 
(V/3 x Vs) 



dY a dY s 



ds 



da 



(Vs x Va) 



= A(a,(9)(Va x V/3) 



(26) 



where A in Eq. (|26|) is constant on field lines (see Eq. ()22[li and depends 
therefore on a and /3 only. This results in the following set of partial differ- 
ential equations 



= \(a,0) (27) 
= (28) 

= 0. (29) 
Vs = |B| lead to 

(30) 



0= B-(R-Y) 
= (Va x V/3) 

+ (R S - F s )Vs] 
= B • Vs(R s - Y s ) 

R s = Y s 





dY a 


da 


d(3 


dY a 


gys 


ds 


da 


dY s 




d(3 


ds 


(Va x 


V/3) 


■[{R a 


~Y a ) 



(R - Y fj )V(3 



(31) 

\B\(R S -Y S ) (32) 

(33) 

due to R — i?°Va + i? /3 V/3 + i? s Vs. The integration procedure is similar 
to that done in Hesse & Schindler (1988). The solution of the differential 
Eqs. 1|28JI and 1)29(1 can be formally expressed by 

d /? s 

■ds + X a (a,f3) (34) 



Y a = 



da 
dR s 
~d[3 



(35) 



so that we get a 'parametric' dependence of the vector Y upon the nonide- 
alness R and therefore a dependence of the new velocity field w upon the 
nonidealness R, in contrast to the formulation in Hesse & Schindler (1988). 
Setting the ansatz Eqs. and into Eq. l(T7|> . we get 
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The formal representation of x a an d X reads 

X a :=J\ a d[3 and X fi ■= J (A + A Q ) da , (37) 

where A" is a function of a and (3 only. From Eas. (|34ll and (|35[) it can 
be seen, that Y a and Y 13 cannot vanish along a field line, as the second 
terms in both equations depend on a and (3 only. If R is strictly localized, 
both components cannot vanish in the direction along the field line, passing 
through the nonideal region. These components are constant on field lines, 
and therefore there will be an additional velocity component to the flux 
velocity, B x Y, which is perpendicular to the magnetic field. But if these 
field lines are not passing through the nonideal region, one can see, that A 
vanishes outside the flux tube. In this case, the field Y could be written as 
a gradient. This enables us to find a common field line velocity for this field 
lines. 

For all field lines passing through the localized (around ao,f3o,so) non- 
ideal region and extending to infinity, the additional velocity component 
B x Y will not vanish, so that the flux conserving velocity field is discon- 
tinuous. 



4 Ohm's law in a partially ionized plasma 
heliosphere 



application to the 



In front of the heliosphere there is a wall of neutral gas, called the hydro- 
gen wall. This leads to the demand, that for describing plasma dynamics in 
this region, it is important not to use ideal Ohm's law, but a more compli- 
cated form of Ohm's law, found by Cowling (1976) and Kulikovskii (1962). 
This was also used by Baranov & Fahr (2003a;2003b), who discussed its 
importance for heliospheric plasma dynamics. 

We will now show, that with Eq. (jf 7(1 it is possible for every point in space 
and at any time to find a flux conserving velocity field because E • B = 
everywhere. 

With the following shape of the nonidealness R 



(I - q ) 2 c 2 
K~i, 



R 



o 



VP x B 



I 

47T 



B x ((V x B) x B) 



Ohm's law can be written as 

E + v x B 
Then Eq. (|39|l can be rewritten as 



R 



E- 



(I - afc 2 



1 + a 



((V x B) x B) 



x B = 0. 



(38) 



(39) 



(40) 



where the term in brackets in front of the cross product with B is the flux 
transporting velocity. But we see, that the velocity depends explicitly on 
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the structure of the (electro-)magnetic field, so it is not clear in advance, 
which should be the precise velocity, and it turns out, that it is not neces- 
sarily a weighted or certain species velocity. This makes it difficult to talk 
about a typical velocity field, in which the magnetic flux is frozen in (see 
e.g. discussion in Hornig & Schindlcr, 1996). So no distinct species velocity 
can be determined as flux conserving velocity. The detection of magnetic 
reconnection in such a plasma is only possible if the physical parameters 
allow to determine a general ideality: there must exist a flux transporting 
velocity almost everywhere. 

5 Discussion and conclusions 

It may be useful finding velocity fields in nonideal environments, in which 
at least the magnetic field lines are frozen in, but maybe not the magnetic 
flux, so that in some regions magnetic flux could be annihilated, diffuse 
or may be even created by special flows (dynamo). The aim is to identify 
velocity fields in a plasma/medium, which includes different ions, neutrals 
and maybe dust, for future work. At this point of the discussion we thus 
can only conclude with the statement that the freezing-in velocity field is a 
very complicated nonlocally determined field depending on many nonlocal 
properties of the field configuration and the differential neutral gas flow 
relative to the plasma flow. 
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